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Abstract. Recent observations of the diffuse Galactic 
gamma-ray glow at 1.809 MeV, attributed to the radioac- 
tive decay of 26 Al, point towards a massive star origin 
of this radioactive isotope. Wolf Rayet stars and core- 
collapse supernovae appear to dominate the production of 
this isotope. Massive stars are commonly located in clus- 
ters and OB associations, regions of recent star formation. 
We thus discuss the temporal evolution of 26 Al, and 60 Fe 
within evolving OB associations. The goal of this study is 
to utilize the associated gamma-ray lines as a diagnostic 
tool for the study of correlated star formation, and also 
to provide more detailed models for the interpretation of 
data obtained with COMPTEL on the Compton Obser- 
vatory. We investigate the effects of possible aluminum 
yield enhancements, predicted for some massive close bi- 
nary systems. In addition to the ejection of chemically 
processed matter, massive stars also drive strong stellar 
winds and emit large fluxes of ionizing radiation. This 
energy and radiation input into the interstellar medium 
(ISM) is crucial for the dynamical evolution of the gas and 
subsequent star forming activity in the galactic disk. We 
discuss population synthesis models for a variety of star 
formation histories, and compare the predicted gamma- 
ray line lightcurves to COMPTEL measurements in the 
Cygnus region. Radioactive tracers such as 26 Al and 60 Fe 
provide a unique gamma-ray tracer of Galactic star for- 
mation activity, complementary to other methods using 
spectral information in the radio, IR or optical bands. 

Key words: stars: abundances, early-type, associations - 
ISM:abundances - 7-rays: observations 



1. Introduction 

For almost a decade the COMPTEL instrument aboard 



many different environments. One of the primary issues in 
the interpretation of this map is the identification of the 
dominant source of this isotope. Several studies of corre- 
lations between the 1.8 MeV map and diffuse maps taken 
in low energy bands suggest that the gamma- ray patterns 
follow closely those of the massive star population (Prant- 
zos & Diehl I1996L Diehl & Timmes |1998l Knodlseder et 



al. 1999). In particular, Knodlseder et al. (1999) recently 



demonstrated that the 1.809 MeV emission profile corre- 
lates best with electron bremsstrahlung emission in the 
microwave regime, which traces the free electrons which 
in turn are produced by the ionizing radiation from mas- 
sive stars. The gamma-ray line from 26 Al thus appears 
to directly trace the presence of massive stars, so that 
Wolf Rayet stars and core-collapse supernovae emerge as 
the m ost promising source candidates (Prantzos & Diehl 
1996Q . With a typical stellar yield of 10~ 4 M Q of 26 Al in- 
jected into the ISM per massive star (during the WR phase 
and in the supernova event) and a million year lifetime, 
one expects a steady state amount of ~ 1 M Q of radioac- 
tive aluminum in the ISM at any time. Spread through- 
out the disk of the Galaxy this amount can explain the 
observed diffuse 1.8 MeV line flux. 

In addition to the nucleosynthesis of radioactive 
matter, which traces the evolution over the past 
few million years, massive stars dynamically shape 
the interstellar medium on a similar time-scale (e.g., 
Lozinskaya 1992). Due to their high mass loss rates 



(Barlow et al. 1981, ie Jager 1988) and large wind ve- 
locities (Cassinelli & Lamers 1987), massive stars impart 
a large amount of momentum and kinetic energy to the 
surrounding medium (e.g., van der Hucht et al. 1987; 
Leitherer et al. 1992). Due to their large surface temper- 



the Compton Observatory (Schoenfelder et al. 1993) has (e.g., Panagia 1973, Vacca et al 



mapped the diffuse Galactic emission in the 1.809 MeV 



aturcs, massive stars also emit a large fraction of their 
radiative luminosity in the wavelength regime below 91.2 
nm, causing photoionization of the surrounding medium 

|1996| ). The subsequent 



line ( Diehl e t al. 1995, Oberlack et al. 1996, and Pliischke 



supernova explosions contribute additional energy (typi- 
cally 10 51 erg (Jones et al. 1998, and references therein)), 



et al. |1999a|) . This gamma-ray line is believed to be due to 
the decay of radioactive 26 Al, which can be produced in 



although much of thi s energy might be radiatively lost 
(Thornton et al. 1998 ), as well as fresh products of stellar 
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nucleosynthesis. 

In a population synthesis approach we compute the 
light-curves of 26 Al (and also 60 Fe) together with the 
mechanical and extreme ultra-violet luminosities as a 
function of time for different star formation histories. 
26 Al and 60 Fe yields are taken from recent WR models 



(Meynet et al. 1997) and supernova simulations (Woosley 
& Weaver |1995| ; Woosley et al. |1995| ). We also study 
the potentially important role of yield enhancements 
for stars that a re me mbers of particular binary systems 
(Langer et al. 



1998|) 



For an assumed star formation 
history the corresponding gamma-ray line flux then 
provides a unique diagnostic tool to study star forming 
regions, complementing other tracers such as IR, UV, 
or H Q emission. We use the gamma-ray light curves to 
constrain the star formation history of the the Cygnus 
region, which is one of the brightest isolated features in 
COMPTEL's 1.809 MeV map. 

The paper is organized as follows. In section |^ we re- 
view the COMPTEL 1.8 MeV results. In section | we 
construct an OB association model based on three major 
aspects; nucleosynthesis of radioactive matter, injection of 
kinetic energy, and emission of extreme ultra violet radia- 
tion. We contrast results from a starburst with those from 
a continuous star formation history. A comparison with 
recent COMPTEL observations of the Cygnus region is 
presented in section |]. We present our concludions in sec- 
tion n 

2. The 1.8 MeV Sky and the Cygnus Region 

The early COMPTEL 1.809 MeV images (Diehl et al. 



1995| ) immediately led to a lively discussion about the 
dominant sources of galactic 26 Al. Radioactive 26 Al is pro- 
duced during hydrostatic H-burn ing in the cores and shells 
of massive stars (Meynet et al. 1997 ), in hydrostatic H- 
shell-burning, in the so-called hot bottom b urnin g (HBB) 
in massive AGB stars (e.g., Bazan et al. 1993 ), during 
hydrostatic Ne-shell-burning in pre-supernova stars, in ex- 
plosive H- and Ne-burning in novae (e.g., Jose et al. 1999|) , 
and core-collapse supernovae (Woosley & Weaver |1995; 



Woosley et al. 1995| ; Thielemann, Nomoto & Hashimoto 



1996; Woosley & Heger 1999). The contributions to alu- 



minum in the ISM from AGB stars, novae, and super- 
novae should lead to different angular patterns of the 1.809 
MeV flux, thus in principal allowing a quantitative reso- 
lution of the flux into the various source populations (see 
Clayton & Leising 1987| , and Prantzos & Diehl |1996| for 



a discussion of this method). The COMPTEL Team re- 
cently published an updated map, based on 8 years of data 



(Phischke et al. 1999a), which confirms the previously re- 
ported characteristics of the galactic 1.8 MeV emission; a 
strong, extended galactic ridge, concentrated towards the 
inner galaxy, a peculiar emission feature in the Cygnus 
region (see Fig. [l]), and a low- intensity ridge extending to- 




Fig. 1. Maximum Entropy image of the 1.8 MeV gamma- 
ray line emission due to the decay of 26 Al in the Cygnus 
region, based on 8 years of COMPTEL observations. The 
contours and colors give the measured 7-ray intensity in 
the line. Over-plotted are known Wolf Rayet stars (dots) , 
supernova remnants (stars) and OB associations (circles), 
wards Carina and Vela. The emerging picture is that mas- 
sive stars and their subsequent supernovae appear to be 
the most promising source candidat es (P rantzos & Diehl 



1996 ). Recently, Knodlseder et al. ( 199S ) carried out an 
extensive comparison between the 1.8 MeV map and a 
variety of other tracers of star formation activity such as 
CO emission, dust maps, and secondary indicators of ion- 
ization (e.g. electron free- free emission). These compar- 
isons strongly support a correlation between massive stars 
and 1.8 MeV emission. Early-type stars are predominantly 
found in regions of recent star formation, so that we can 
test the global correlation with a detailed analysis of the 
Cygnus star forming region. Fig. [I] is an enlarged image 
of this region drawn from the COMPTEL 1.8 MeV sur- 
vey. Due to their high mass loss rates Wolf Rayet stars 
and core-collaps supernovae are primary candidates for 
26 Al sources. In Fig. [I] Wolf Rayet stars and supernova 
remnants are marked by dots and stars, respectively. The 
encircled regions mark identified OB associations in the 
Cygnus complex. The Cygnus region includes at least 21 
Wolf Rayet stars and 17 supernova remnants, as well as 8 
OB associations of different richness class and age. 

3. Modelling OB Associations 

The Wolf Rayet phase and the core-collapse supernova ex- 
plosion are specific episodes during the evolution of mas- 
sive stars. Therefore groups and associations of early-type 
stars are key targets for 26 Al abundance studies. Massive 
stars have a significant dynamic impact on the surround- 
ing interstellar medium, especially when their winds and 
explosions occur in near- simultaneity Star formation is 
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known to occur in groups with a hierarchical size distri- 
bution (e.g. Elmegreen, D.M. & Salzer, J.J. |1999|, Efrc- 



mov, Y.N. & Elmegreen, B.G. 1998), ranging from small 
groups to spiral arm fragments. The mini star-burst lead- 
ing to group or cluster formation might lead to a coeval 
population of ZAMS stars, but it might also be spread 
over some significant span in time. But the subsequent 
stellar wind and supernova phases are still significantly 
localized (in space and time), so that the energy input 
into the ISM from multiple explosions can lead to the for- 
mation of super-supernova remnants, c ommo nly referred 
to as superbubbles (e.g. Tomisaka et al. 1981 ). We there- 
fore study the properties characterizing the interactions 
of groups of massive stars with their surrounding medium 
with a self-consistent population synthesis approach. 
Our aim is to establish a numerical model, based on avail- 



able st ellar e voluti on models (Maeder et al. 1994 ; Woosley 
et al. 1993 , 1995 ), of the temporal evolution of a star 



forming region, with emphasis on gamma-ray line pro- 
ducing isotopes. To achieve a simple and transparent de- 
scription of the investigated parameters, the model uses a 
semi-analytic approach. We utilize the detailed models on 
the basis of simple fit- functions; e.g. power-law descrip- 
tion for the stellar lifetime T star as function of the initial 
mass Mi. All our simulations start with solar composi- 
tion, but incorporate subsequent chemical evolution. The 
inclusion of core-collapse supernovae is based on parame- 



ter fits to recent calculations by Woosley & Weaver (1995) 



and Woosley, Langer & Weaver (1995) 



Stellar groups are fundamentally characterized by their 
initial mass function $ (IMF) and star formation rate S 
(SFR). We adopt the common assumption that the IMF is 
time- independent and that the SFR is mass- independent. 
The model employs a lower mass limit of 8 Mq, thus re- 
stricting the simulation to stars that end their life in a 
core-collapse supernova. We adopt a single power-law IMF 
(eq. 0) with exponent T = 1.35, which is relevant for high 
mass stars considered here (Kroupa 2000 ; Scalo 1998 ). The 
normalisation constant a was chosen such as to normalize 
the IMF to unity for the mass range under investigation, 
so ^(Mi) gives the probability of existence of a star with 
initial mass Mi. 



*(M 4 ) = a ■ M" (1+r) 



(1) 



As star formation histories we consider either an instan- 
taneous burst or a continuous, constant star formation 
rate with an adjustabe duration. We verified the simu- 
lation through consistency checks against model calcula- 



tions by Leitherer et al. (1992) and results from the STAR- 
BURST99 code (Leitherer et al. |1999| ). Stellar evolution- 
ary phases of particular relevance in our model are the 
main sequence/post main sequence phase, the Wolf Rayet 
phase, and the terminal core-collapse supernova. In the 
case of 26 Al synthesis we furthermore include the effects of 
the presence of massive close binary systems (MCB). For 
the calculation of the stellar wind properties we specifi- 
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Fig. 2. Schematic overview of the major components of 
our OB association model. 

cally incorporate the BSG/RSG and LBV phases. Fig. 
provides a schematic overview of the model ingredients. 
We now describe in greater detail the different processes 
included in our model, beginning with the ejection of ra- 
dioactive isotopes with intermediate lifetimes. Then we 
discuss the description of matter ejection and kinetic en- 
ergy flux, and present our estimates of the ionizing Lyman 
continuum emission. 



3.1. Synthesis of the Radio-Isotopes Al and Fe 

In addition to 26 Al with a mean lifetime of r = 1.04 Myr 
and an associated 7-ray line at E = 1.809 MeV, we also 
consider the important radioactivity 60 Fe, with a lifetime 
of 2.07 Myr and comparable yields in core collapse super- 
novae. 60 Fe gives raise to two 7-ray lines in the COMPTEL 
energy regime at 1.173 MeV and 1.332 MeV. Whereas 
26 Al can be synthesized in hydrostatic as well as explo- 
sive scenarios, explosive nucleosynthesis is the dominant 
production mechanism for 60 F e (Ti mmes et al. 1995 and 
refs therein; Diehl & Timmes 1998). As was emphasized 
by Timmes et al. 26 Al and 60 Fe are co-synthesized in the 
same shells of the collapsing massive star, so the explosive 
yields of these nuclides are expected to be well correlated. 
26 Al is produced hydrostatically via proton capture on 
25 Mg. In this scenario, 26 Al is predominately produced 
during hydrostatic hydrogen burning in cores and shells 
of massive stars. To be seen by its 7-ray emission, 26 Al 
has to be expelled into the ISM by the star. Thus, effec- 
tive mixing as well as significant mass loss due to stellar 
winds are further prerequisites for significant 7-ray emis- 
sion from hydrostatic 26 Al nucleosynthesis. Due to their 
high mass loss rates, Wolf Rayet stars are expected to be 
significant contributors to the steady state abundance of 
26 Al's in the ISM, and some WR stars are near enough to 
be potentially detectable as point sources. Fig. || shows the 
ejection rate (mass per unit time) from two representative 
stars as a function of time. The transition phase of early 
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Fig. 3. 26 Al emission (in solar masses per year) as a func- 
tion of time for a 25 M Q (upper panel) and a 120 M Q 
(lower panel) star. 



type stars (Of) towards the post main sequence phases 
and the early Wolf Rayet phase turn out to be the dom- 
inant ejection pha ses fo r 26 Al. 26 Al ejection from AGB 
stars (Bazan et al. 1993 ) and nova explosions (Jose et al. 
19991 ) is neglected in this study for two reasons. First, the 



COMPTEL 1.8 MeV all sky map correlates best with trac- 
ers of the young population of massive stars, as discussed 
in section ^. Second, the timescales involved in AGB and 
nova phenomena greatly exceed the typical lifetime of OB 
stars and their associations. 

Thus, in modelling 26 Al and 60 Fe nucleosynthesis in OB 
associations we focus on contributions from hydrostatic 
nucleosynthesis of massive stars and explosive synthesis 
in core-collapse supernovae. Whereas in the case of SNe 
the ejection of nucleosynthesis products is a singular event 
in time (eq. ||), the release of hydrostatically produced 
species by stellar winds has to be treated with a time- 
dependent model. This is of particular importance for 
stars for which the duration of 26 Al ejection is of the 
same order, or even much longer, than the lifetime of 
26 Al. Fig. H shows two representative temporal profiles 
for stars with initial masses of 25 and 120 Mq, respec- 
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Fig. 4. 26 Al and 60 Fe lifetime-integrated yields as a func- 
tion of the stellar mass on the ZAMS. 



tively. These profiles are derived from detailed evolution- 
ary models of Maeder & Meynet ( 1994[ ) and Meynet et al. 
(1997). To restore the analytic character of the model, we 
apply a Gaussian approximation (eq. 0) for these profiles. 



The error of this approximation is less than 10%. The pa- 
rameters of the Gaussian (area, position and width) are 
computed on the bases of interpolating polynomials (e.g. 
area(Mi) — P^(Mi)). Contributions from explosive nu- 
cleosynthesis are calculated on the bases of core-collapse 
supernovae simulations (Woosley & Weaver 1995| ; Woosley 
& Heger 1999; Woosley, Langer & Weaver 1995] ). A direct 
comparison of explosive yields given in Woosley & Weaver 



(1995) with results from more recent models including ro- 



tation and mass loss (Woosley & Heger 1999) gives good 
confidence in the older results. Thus we use the yields from 
the 1995 paper of Woosley & Weaver, covering a wider 
mass range, to model type II supernovae. Type Ib/c su- 
pernovae from more massive stars, showing a Wolf-Rayet 
phase, are more d elicate to incorporate. Woosley, Langer 
& Weaver ( |l995|) have calculated supernova models for 
mass lossing He stars. Supernova yields presented by this 
group have to be care fully associated (via the CO core 
mass (cf. Maeder 1992| )) to evolutionary tracks of model 
stars from Maeder et al. ( |1994 ). Unfortunately, these su- 
pernova models only cover CO core masses significantly 
lower than those resulting from mass lossing, massive stars 
(Mzams > 40 Mq). On the other hand the CO core 
masses from the massive stars show only minor variations 
as function of the initial mass and lie in the range of a 
25 Mq star. Thus we decided to model the radio-isotope 
yields from type Ib/c supernovae by a constant value. Fig. 
[I| displays the lifetime- integrated 26 Al / 60 Fe yields as func- 
tion of the initial mass. In the case of 26 Al the integrated 
yield for stars more massive than 40 Mq is strongly dom- 
inated by the wind ejected 26 Al. In the case of 60 Fe the 
wind ejected yield is negligible. 

To calculate the amount, M x (t), of radio-nuclides x we 
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Fig. 5. 26 Al lightcurves from an OB association at a 
distance of 1 kpc in the case of a burst-like star forma- 
tion converting 10 4 Mq into stars (upper panel) and a 20 
Myr episode of constant star formation rate that forms 
the same mass of stars (lower panel). 



solve equation with the approximations given in equa- 
tions m and H (see text above). 



M x (t) = 



dtZ(i) / dMi$(Mi) 

3 JM min 



diy x {Mi,t-t) -exp 



t-t 



y h x vd (t) 



er(Afj) v 27r 



exp 



{t-t Peak {M t )f 



2a 2 {Mi) 

»T (t) = £f v (Mi) ■ S((t - t) - w (Mi)) 



(2) 

(3) 
(4) 



We calculate M 2 e(t) as well as M 60 (t) for various sets of 
initial mass function and star formation rate histories. 
Fig. H and ^| show the resulting 26 Al and 60 Fe lightcurves 
for an instantaneous starburst and a case with continuous 
star formation, respectively. In both cases the total mass 
of stars created was assumed to be 10 4 Mq, with an IMF 
exponent of T = 1.35. The duration of the star formation 
activity in the contiuous case was 20 Myr. To convert 
the resulting decay rate to gamma-ray line fluxes we 
assumed a nominal distance to the star forming region 
(OB association) of 1 kpc. Lightcurves are calculated 
using four different upper mass limits in the IMF; 60, 
80, 100, and 120 Mq. A comparison of the lightcurves 
shows a strong dependence of shape on the assumed 
star formation history. In particular, the double-peaked 
shape of the burst lightcurve is completely washed out 



Fig. 6. Same as Fig. |^ for the sum of two high energy 
7-ray lines of 60 Fe at 1.173 and 1.332 MeV. 



in the continuous scenario. Furthermore, peak fluxes are 
reduced by a significant factor up to 5. 
Using COMPTEL's line sensitivity of 10~ 5 phcm~ 2 s _1 , 
we derive upper limits for the detection distances to such 
regions of star formation: 2.4 kpc at maximum of 26 Al 
for the starburst scenario, and 1.3 kpc in the continuous 
case, whereas for 60 Fe the corresponding values are 1.5 
kpc and 1 kpc. 

Fig. shows the time-evolution of the 60 Fe to 26 Al 
flux ratio. After an initial increase the ratio remains 
constant at ~ 10-20%. The increase at late times is due 
to the vanishing hydrostatic component of 26 Al emission 
and the longer lifetime of 60 Fe. In the continuous case 
associations with an smaller upper mass limit show an 
increased ratio, due to the suppression of additional 26 Al 
ejection from very massive Wolf Rayet stars. Only during 
the late phases this ratio increases due to the lifetime 
difference. 

Stars are often observed to form binary systems (Mason 
et al. 1998. Depending on the relative orbital parameters, 
mass-transfer reactions via Roche-Lobe overflow (Kfiz 
1970) may accure. These binary systems are labeled 
'close' systems. In massive, close binary systems mass 
transfer is expected to significantly alter the evolution of 
the individual stars. These modifications may also induce 
significant variations in the nucleosynthesis and the sub- 
sequent enrichment of the surrounding ISM compared to 



single stars. Langer et al. (1998) studied 26 Al nucleosyn- 



thesis in particular massive close binary systems, based 
on evolutionary models that include mass transfer from 
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Fig. 7. Flux ratio of 60 Fe to 26 Al as function of time for a 
starburst scenario (top) and a continuous star formation 
(bottom). The 60 Fe 7-ray flux is estimated as in Fig. |(| 



Fig. 8. Ejected 26 Al mass as function of time for OB as- 
sociations with MCB component (top: starburst scenario/ 
bottom: continuous SF) 



the primary star. These simulations are still somewhat 
limited in scope, due to the neglect of stellar winds 
during the evolution of the massive stars. Furthermore, 
the calculations were only carried out for quarter solar 
metallicity. Still, these preliminary investigations strongly 
suggest that the yield of 26 Al can be greatly enhanced 
during the SN explosion of the secondary, perhaps by as 
much as a factor 10 3 . The origin of such a large enhance- 
ment is attributed to a 'positive interference' of several 
effects caused by mass-transfer during the main-sequence 
evolution of the orbiting stars. The models presented by 
Langer et al. assume case A or case AB mass-transfer, 
which labels mass-transfer via Roche-Lobe overflow 
during the core hydrogen burning phase of the primary 
star (Kippenhahn & Weigert 1967 ). Due to the loading 
of fresh material on to the outer layers of the secondary 
star, the secondary star is rejuvenated. Depending on the 
details of the mixing processes the opacity of the outer 
layers of the secondary star is reduced, which causes the 
mass loss rate to drop significantly. Stars modified in this 
manner exhibit altered evolutionary tracks (i.e. no WR 
phase) and end their lives in a supernova explosion as a 
BSG-likc star (similar to SN1987A). Due to the smaller 
evolutionary gap between the main-sequence and BSG 
phase, the time-delay between 26 Al production and its 
explosive ejection is drastically reduced. This can lead 
to a huge enhancement of the yield. The adoption of 
the modified MCB yields to the framework of our model 
OB associations, which must rely on solar metallicity 
models, consequently introduces significant uncertainties. 



For stars of solar metallicity the neglect of stellar winds 
during the evolution of the individual stars is much more 
critical than for the case for which Langer et al. worked 
out their model. We therefore decided on a conservative 
approach, limiting 26 Al yield enhancements to factors 
less than 10 2 . MCB systems are included in the OB 
association model with a typical frequency of ~_5%. 
This number is based on results of Mason et al. ( 1998 ), 



who finds the total binary frequency in clusters and 
associations in the solar n eighbo rhood to be at least 59%, 
while Maeder & Meynet fll994j ) argued that 8% of such 
binary systems are massive star systems that undergo 
a mass transfer phase. We also assumed that all MCB 
systems show the same enhancement effects. Furthermore, 
we introduced an upper mass limit for stars in MCB 
systems in the range 40 to 60 Mq. Fig. || shows the 
expelled 26 Al mass, M^eit), as a function of time. Except 
for the additional MCB component all parameters of the 
simulations are identical to those used in the previous 
cases. Given the large enhancement factors it comes to 
no surprise that the M26 time profiles are now dominated 
by the population of MCBs. The SN peak is much more 
pronounced than in the case of associations consisting 
exclusively of single stars. In addition, a prolonged 
evolutionary time in mass-transfer systems causes a time 
shift that alters the shape of the aluminum profile. The 
gray shaded area marks the uncertainty due to lifetime 
variations. Fig. ^ displays the 26 Al overproduction due 
to MCBs, in comparison to a pure single star association 
that underwent a starburst (keeping the total number 
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Fig. 9. The gray lines give the 26 Al timcprofiles for an 
upper mass limit of 120 Mq for the two different popu- 
lations. The black lines show the overproduction of 26 Al 
due to enhanced yields in an association containing MCB 
systems. 



of stars identical in both cases). Due to uncertainties 
in the MCB implementation, our results for an OB 
association with an MCB population should be considerd 
preliminary. Further studies of nucleosynthesis in MCB 
systems that include stellar winds are required to proceed 
in this investigation. 

3.2. Kinetic Energy 

Massive stars and supernovae are strong sources of kinetic 
energy, which significantly affects the dynamics of the 
surrounding interstellar medium. During the Wolf Rayet 
phase massive stars experience serious mass loss rates 
of order lO^MQ yr" 1 flBarlow et al. 1981| , |Doom 1988 , 



Nugis et al. 1998). This matter is expelled into the sur- 
rounding me dium with velocities of ~ 10 3 kms _1 (e.g., 
Prinja et al. |l990| ). The corresponding mechanical lumi- 



nosities therefore reach values of order 10 37 erg s _1 for indi- 
vidual Wolf Rayet stars. The contributions from OB stars 
during the main sequence and red/blue (RSG/BSG) su- 
pergiant phase are considerably smaller, due to lower mass 
loss rates and/or lower wind velocities. The contribution 
from LBV (luminous blue variables) stars are small as 
well, due to the short duration of this evolutionary phase. 
Still, we included these phases in our model. In OB as- 
sociations continuous energy injection from stellar mass 
loss coincides with punctuated energy injection from su- 
pernovae. 

Supernovae typically release 10 51 erg of energy into the 
surrounding medium (e.g., Jones et al. 1998). With a su- 
pernova rate in the range of 10 -6 to 10 -4 per year one thus 
finds an average mechanical luminositiy from SNe that is 
comparable to that from stellar winds of massive stars. 
We adopted mass loss rates given in (Maeder & Meynet 



1994) and wind velocities calculated from a fit function 
provided by Leitherer et al. ( |1992| ). For the supernova 
contribution we assume the standard value of 10 51 erg per 
event for all events. 
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Fig. 10. Timeprofiles of the mechanical luminosity for an 
OB association in a starburst as well as continuous star 
formation scenario. 

We calculate the timeprofile of mechanical luminosity via 
equation where l w is the mechanical wind luminosity of 
a single star as function of initial mass and age. For two 
limiting cases (burst vs. continuous star formation) Fig. 
|To| displays the resulting timeprofiles. 

dN stars (t) 



L w (t) 



■ 10 51 erg 



+ 



dtY>(i) 



In 



dMi<Z>(Mi) 



dtl w (Mi,t- i) 



(5) 
(6) 



l w (M,t) - --M(M,t)-v 2 w (M,t) 

In the starburst scenario the wind and the supernova 
phase can easily be distinguished. During the early evo- 
lution of a burst-formed OB association, stellar winds are 
the only sources of kinetic energy. Whereas in the later 
phases supernova events take over and balance the me- 
chanical luminosity near 10 39 ergs _1 . In the case of con- 
tinuous star formation the mechanical luminosity is al- 
most constant over a time period of up to 15 Myr. Then 
it slightly decreases due to a decrease in the stellar wind 
contribution. 

In addition to the mechanical luminosities we also calcu- 
late the total mass loss rate. The mass liberated during 
the supernova events is determined by assuming that each 
SN leaves behind a compact remnant with a mass of 1.4 
Mq. Fig. [fl] shows the mass loss rate (integrated over the 
population) as a function of the age of the association for 
the two limiting cases of the star formation history. The 
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Fig. 11. Mass loss rate of an OB association for a star- 
burst and continuous star formation, respectively. 




Fig. 12. Timepronles of photon number intensity for pho- 
tons with wavelengths shorter than 91.2 nm. 



simulations are based on an identical parameter set to 
that described in the nucleosynthesis section, but without 
inclusion of MCB systems. A comparison between stellar 
wind and supernova mass expulsion shows that the stellar 
wind is the dominant source of matter blown into the ISM. 
However, one should keep in mind that ejecta processed 
in a supernova are more metal enriched than the matter 
ejected by stellar winds. Thus, supernovae will always in- 
crease the metallicity of the surrounding medium. 

3.3. Lyman Continuum Emission 

Besides the mechanical energy, massive stars also release 
a large amount of energy in the form of radiation. Due to 
their high effective tem perat ures (~ 10 4 K) O stars (e.g. 
see Maeder & Meynet 1994 ) emit a significant portion 
of their electromagnetic radiation in the wavelength 
regime below A = 91.2 nm where hydrogen, the by far 
most abundant element in the interstellar medium, is 
photoionized. We included this ionizing radiation in 
our model by fitting the stell ar Lyman photon fluxes 
estimated by Vacca et al. 1996 The adopted fit-function 
(eq. 0) is accurate within 5% for the whole mass range. 



Q 0/1 (Mi) = exp ai 



02_ 

Mi 



-r=r • 10 



49, 



(7) 



One major drawback of this approach is the loss of spec- 
tral information. We only compute the number intensity 
of photons capable of ionizing hydrogen (Qq) and helium 
(Qi). Fig. |lj shows the logarithmic temporal profiles of 



these integrated EUV-photon fluxes. In the starburst sce- 
nario the ionizing flux decreases very quickly, whereas in 
the case of continuous star formation a plateau phase is 
quickly reached, and lasts essentially for the whole dura- 
tion of the star formation activity. 

4. 26 Al FROM CYGNUS OB ASSOCIATIONS 

As mentioned in section 0, the Cygnus region provides 
a suitable site for a detailed comparison between model 
predictions and data. Besides the many likely 26 Al point 
sources in this region (Wolf Rayet stars and supernova 
remnants) Cygnus contains at least 8 OB ssociations. 
We use our OB association model to interpret the recent 
COMPTEL 1.8 MeV results from this field. First results 
from such a comparison were presented by Pliischke et al. 



(19991). In a sub-sequent paper we will summarize our 
analysis of the Cygnus region in detail, for now we only 
want to summarize shortly the main results from analysis 



given in Pliischke et al. (1999b). 

On the basis of age and richness estimates by Bochkarev & 
Sitnik ( |1985| ) and Massey et al. ( |1995D we normalized the 
IMF by adopting a single power-law with mean exponent 
r = 1.1 (Massey et al. 1995) (see eq. |l|). We then compute 
26 Al lightcurves for starburst scenarios as well as Gaus- 
sian star formation histories with a maximum duration of 
15% of the estimated age of the association. The result- 
ing gamma-ray line fluxes have significant uncertainties 
due to possible variations in the IMF and especially the 
uncertain MCB contribution. While the model can repro- 
duce the observed flux from the Cygnus region, there is a 
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Fig. 13. Model/Measurement comparison for the Cygnus 
region. 



possible variation of a factor ~ 5, which renders the inter- 
pretation somewhat inconclusive. Additional constraints 
(x-ray flux, dynamic state of the gas in the star form- 
ing region, etc) have to be used to improve the model and 
thus refine the constraints derived from it. Fig. |l3| shows a 
coarse model-measurement comparison. The gray shaded 
image in the upper panel is a maximum entropy image 
of a model calculation including MCBs whereas the con- 
tour lines mark the COMPTEL maximum entropy result. 
The lower panel shows a direct comparison of the inten- 



sity profiles, taken between the dashed lines in the map. 
The distributions along the galactic plane are comparable, 
whereas the reconstructed fluxes are somewhat too low in 
this example. In future studies we shall incorporate ad- 
ditional diagnostics of the model to refine the analysis of 
the Cygnus region. 

5. SUMMARY & CONCLUSION 

We discussed a simple population synthesis model for de- 
termining the evolution of the abundance of gamma-ray 
line producing radiactivity in OB associations. We also 
traced the stellar output of kinetic/radiative energy and 
ionizing luminosity, which are essential for studies of the 
impact of OB associations on the surrounding ISM. A cor- 
nerstone of our model is the prediction of emission prop- 
erties of radioactive isotopes with lifetimes comparable to 
the main-sequence lifetime of the massive stars that dom- 
inate the dynamical evolution of OB association. In par- 
ticluar, we focus attention on 26 Al and 60 Fe, because the 
gamma-ray lines from these two species provide a poten- 
tially powerful tracer of the star formation history in the 
association. We advocate the use of gamma-ray line obser- 
vations as a complementary diagnostic tool. Beside the nu- 
cleosynthesis of radioactivities we also model the injection 
of kinetic energy and extreme ultra-violet radiation. For 
our future work on the evolution of OB association we plan 
to use the mechanical luminosity (the central driver) to 
combine gamma-ray line tracers with a dynamical model 
for the evolution of a super-bubble, driven by the winds 
from massive stars and multiple SN events. Based on the 
results presented here, we draw the following conclusions: 

1. Massive stars and their supernovae will give raise to a 
dynamic evolution with significant temporal variations 
of the abundances of 26 Al and 60 Fe. The resulting time 
profile of the associated gamma-ray lines can be used 
to constrain or even determine the star formation his- 
tory of a stellar cluster (OB association). 

2. In steady state the 60 Fe to 26 Al ratio will be of the 
order 10 - 20% (e.g., Timmes et al. 1995), but in an 
evolving OB association, this ratio provides an excel- 
lent constraint on the age and thus dynamical stage of 
the system (see Figure 7). 

3. MCB systems could dominate the shape of the 26 Al 
lightcurve, but the study of the nucleosynthesis in 
these systems is not yet mature enough to draw final 
conclusions. 

4. The mass of injected matter into the local interstellar 
medium due to stellar phenomena is strongly domi- 
nated by stellar winds (see Figure 11). 

5. The short life-times of early- type stars cause a very 
rapid decrease of LeuVi so EUV emission is impor- 
tant in very early evolution phases of OB associations. 
This is also the time when the gamma-ray lines are at 
detectable levels. It is thus important to consistently 
treat tracers based on the emitted ionized radiation 
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and the gamma-ray lines from radioactive aluminum 
and iron. 

6. Detectibility of 60 Fe: In the 1 MeV regime the narrow 
line sensitivity of the INTEGRAL spectrometer SPI 
is expected to be nearly one order of magnitude lower 
than COMPTEL's sensitivity (for a on-axis observa- 
tion of 10 6 s). Based on the presented analysis of the 1.8 
MeV emission of the Cygnus region and the prediction 
of the flux ratios (see Figure one expects Cygnus 
to be visible in the 60 Fe lines using INTEGRAL. 



A comparison of our predictions for 26 Al with recent mea- 
surements in the Cygnus region shows that our model is 
capable of reproducing the measured fluxes. On the other 
hand, the comparison also revealed the need for additional 
observables to better constrain the models. The inclusion 
of enhanced aluminum production in MCBs emphasizes 
the possibility of large uncertainties due to additional free 
parameters. 

The study of star formation in the Galaxy and through- 
out the universe is a fundamental pillar of modern astro- 
physics. We know that the global SFR changes by more 
than an order of magnitude over cosmic time. The Milky 
Way galaxy also underwent significant evolution of its SF 
history (though not congruent to that of the universe as a 
whole), but we do not know the exact shape of the SFR(t) 
function. We do not even know the exact value of the 
galaxy wide SFR at the present time, but is of order of a 
few solar masses per year. We also know that the star for- 
mation activity is not spread uniformly though the disk, 
but instead located in a hierarchy of star forming cen- 
ters (groups, clusters, a ssocia tions, complexes, etc., e.g. 
Elmegreen & Elmegreen 1996 ). The strong correlations in 
space and time of the SF process lead to many dynam- 
ically active cradles of star formation that in turn drive 
the dynamics of the ISM and the Galactic disk as a whole. 
The study of these star forming regions is thus vital to our 
understanding of galaxy evolution on a larger scale. In this 
paper we promote the use of gamma-ray lines from radi- 
active aluminum and iron to study the SF process in these 
cradles, as well as galaxy wide. The gamma-ray map does 
not suffer from extinction, so that we get a good global pic- 
ture of ongoing star formation throughout the Milky Way. 
Although the flux maps do not contain direct information 
of distance, realistic spatial models of the emission suggest 
that the Galaxy indeed converts a few solar masses of gas 



into stars each year (Timmes et al. 1997). Here we have 
focused on the Cygnus region to the demonstrate how the 
gamma-ray line method can be applied to individual star 
forming complexes. The use of radiactive tracers is a rela- 
tively recent addition to our arsenal of probes of SF, one 
that complements other more traditional tools and that 
should be included in any study of recent star formation 
activity in the Galaxy. Of course, the life times of the 
isotopes promoted here limits applications to star forma- 



tion within the past 10 Myrs, unless one considers nearby 
sources. 
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